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ABSTRACT 

The dynamics of relativistically strong electromagnetic (EM) wave prop- 
agation in hot electron-positron plasma is investigated. The possibility of 
finding localized stationary structures of EM waves is explored. It is shown 
that under certain conditions the EM wave forms a stable localized soliton- 
like structures where plasma is completely expelled from the region of EM 
field location. 
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During the last few years a considerable amount of work has been de- 
voted to the analysis of nonlinear electromagnetic (EM) wave propagation in 
electron-positron (e-p) plasmas [1]. Electron- positron pairs are thought to 
be a major constituent of the plasma emanating both from the pulsars and 
from inner region of the accretion disks surrounding the central black holes 
in the active galactic nuclei (AGN) [2] . The process of e-p pair creation oc- 
curs in relativistic plasma at high temperatures, when the temperature of the 
plasma becomes of the order of, or larger than, the rest energy of electrons. 
Such relativistic plasmas have presumably appeared in the early universe 
[3]. Intense relativistic e-p plasmas could also exist in the vicinity of cosmic 
defects hke superconducting cosmic strings [4]. Collective processes in e-p 
plasmas are of considerable interest. The processes of wave self- modulation 
of EM waves and soliton formation have attracted a great deal of attention. 
Stable localized solution may be a potential mechanism for the production 
of micro-pulses in AGN and pulsars [5]. In the early universe stable localized 
EM waves could create inhomogeneities necessary to understand the observed 
structure of the visible universe. 

In the recent paper of Berezhiani and Mahajan [6] the nonlinear propa- 
gation of relativistically strong EM radiation in a hot e-p plasma has been 
considered. It has been shown that e-p plasma supports the propagation of 
nondiffracting and nondispersive EM pulses (hght bullets) with large den- 
sity bunching. However, the authors concentrated in the case of transparent 
plasma and consequently the group velocity of the pulses is close to the ve- 
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locity of light c. In the present paper we consider the propagation of strong 
EM radiation in a hot unmagnetized e-p plasma aiming to find the localized 
stationary soliton-type solutions. 

We start from Maxwell equations to describe the EM wave propagation 
in an e-p plasma, expressing the fields by the vector and scalar potentials, 
i.e. 

1 dA 

E=— — -V(^ , B = VxA (1) 

c at 

where the coulomb gauge V • A = is fulfilled. Accordingly the field equa- 
tions take the form: 

^ - A A + ( Vy.) - AttcJ = (2) 

and 



Acp = -47rp (3) 
Here, p and J are the charge and current densities given by 

P^Y. ^oc-ria, J = H (^ana^a (4) 

a a 

where a indicates the particle species a (= e.p for electrons and positrons, 
respectively); and Ha are the charge and density of the corresponding par- 
ticle a. We consider the case in which the equilibrium state is characterized 
by riQe — riop — Uq, where rioa is the equilibrium density of the particle a. 

Before writing the hydrodynamic equations of relativistic plasma it is 
necessary to define what relativistic means. In fact, we have two types of 
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relativistic regimes in plasma: in a strong EM field the plasma particles may 
obtain relativistic velocities. In space, the EM radiation of objects (nuclei 
of galaxies, radio-galaxies, quasars, pulsars, etc.) may serve as a source of 
such strong fields. The case when the thermal energy of the plasma particles 
is of the order of, or larger than, the energy at rest, it is the other type of 
relativistic regime. In this case the thermal velocities of the particles become 
of the order of the light speed. Certainly, in both cases the decisive role 
belongs to relativistic effects in plasma, but the character of its manifestation 
is different. Both relativistic effects can play an important role in the e-p 
plasma. Let us assume that the velocity distribution of the particles of species 
a is locally a relativistic Maxwellian. Then, according to Ref.[6, 7], the set 
of relativistic hydrodynamic equations of motion can be written as: 



^(moaG'a7aC^) " —^Pa = CaU^E (5) 

at Ur, at 



^APaGa) + — VP„ = e„E+ X B) (6) 

dt Ua C 

The continuity equation for the particle a is 

^ + V(n„u„) = (7) 

ot 

Here pc = ■ja'iTT'Qa^a IS the hydrodynamic momentum. Pa = UaTa/ja is 
the relativistic particle pressure, is the hydrodynamic velocity of the fluid, 
7a = (1 — ul,/ c^)"^^'^ is the relativistic factor, moa and are the particles 
invariant rest mass, and temperature respectively, d^/dt = d/dt + Uq,V 
is the comoving derivative. The role of the particle- mass is now played by 



the quantity Me// = moaGa{za) , where Ga{za) = K^{za) / K2{za) ■ Here 
Kz{za) and ^2(^0) are respectively the McDonald functions of the second and 
third order [z^ — moaC^/Ta) . The effective mass of the particles Mg// de- 
pends on the temperature. For nonrelativistic temperatures (Tq, <^ moaC^) 
Mf,ff = moa(l + 5Ta/2mQaC^) , while for ultrarelativistic high temperatures 
{Ta ^ niQaC^) the effective mass becomes Mg// — ATa/c^ , and the fluid 
inertia is primary provided by random thermal motion of the particles. In 
this case, the rest-mass is negligibly small and the e-p gas behaves like pho- 
tons. Using simple manipulation, from the Eqs. (5) and (6), we obtain the 
adiabatic equation which reads 

exp{-ZcGa(zc,)) = const. (8) 

In the nonrelativistic limit {za ^ 1) , Eq. (8) yields the result for a 
mono atomic ideal gas {ua/ (Tq^^q^^) — const.) , and in the ultrarelativistic 
limit {za <^ 1) we have the adiabatic law for the "photon" gas {ua/i'jaT^) — 
const.). 

We are looking for a localized one-dimensional solution of this system of 
equations for a circularly polarized EM wave, where the vector potential A 
can be expressed as: 

A_L = -(x -|- i'y)A{z, t)exp{—iuot) + c.c (9) 

where A{z,t) is a slowly varying function of t, cuq is mean frequency , x and 
y are the unit vectors. The transverse component of equation of motion (6) 
is immediately integrated to give (for details see Ref. [6]): 
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Pa±Ga = --A^ (10) 

c 

where the constant of integration is set equal to zero, since the particle hy- 
drodynamic momentum is assumed to be zero at infinity where the fields 
vanish. 

Before writing the equations for the longitudinal motion we would like 
to mention that this motion of the plasma is driven by the ponderomotive 
pressure (~ p^j^) of high-frequency EM fields and it does not depend on 
the particle charge sign. In what follows we assume that in equilibrium the 
temperatures of electrons and positrons are equal, i.e. Tq^ — Tgp — Tq. Since 
the effective mass of the electrons and positrons are equal {Gg — Gp — G), 
the radiation pressure gives equal longitudinal momenta to both the electrons 
and positrons (pez = Ppz = Pz) and affects concentration without producing 
the charge separation. Consequently Ue — — n and = 0. It is also 
evident that due to symmetry between electron and positron fluids their 
temperatures remain equal {T^ — Tp = T). 

It is now convenient to introduce the following dimensionless quantities: 



Pa n 1 eA CJe + 

Pa = , n=—, T = -, A = -, r = — r, t = (jj^t (11) 

moe no moeC^ rriQeC^ c 

where = (47re^no/'n2oe)^^^ is the electron Langmuir frequency. 

The longitudinal motion of the plasma is determined entirely by the set 
consisting of the z component of the equation of motion (6), 
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and the "energy" conservation equation (5), 

f d 1 dnT 1 d\A\^ 

where = Pz/l- The relativistic factor 7 does not depend on the "fast" 
time (<^o'^) and can be written as: 



7 = 



(14) 



G2 

Substituting (9) and (10) into Eq.(2), then for a slowly varying amplitude 
of EM wave A{z, t) we obtain the following equation: 



dA d'^A 

'^i^0-g^ + -Q;2+^-^ + '^f^-^ (15) 



where 



_ nGo(To) 

7G(T) ^ ^ 

and A = uJq — 2. For convenience we redefined the electron rest mass in 
Eq.(ll) as moe — > ?77.oeGo(3o). In dimensional units A ~ a;Q — 2^"^, where 

coe = (47re2no/moeG'o(ro))'/'. 

We are looking for the stationary localized solutions (vanishing at infinity) 
described by Eqs. (8), (12)-(16). Assuming that \A\ depends only on the 
special coordinate z, and integrating Eqs. (12) and (13), we get the following 
integral of motion 
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G{T)^ = Go{To) 



(17) 



where 7 = (1 + 1^17^2)1/2 ^ Ymva Eq. (17) we get: 



A 



2\ 1/2 



G = Gn 1 - 




(18) 



It follows from Eq.(18) that the present hydrodynamical theory, which 
describes the nonhnear waves in e-p plasma, is vahd for |74|^„^/Go < 1. 
When the latter is violated, then the electromagnetic waves are overturned 
and they will cause multi-stream motion of the plasma. In such a situation, 
one must resort to kinetic description for studying the nonlinear wave motion. 
This investigation is, however, beyond the scope of the present paper. 

Using Eqs. (17)-(18) we get f = 1 — n. Now we should obtain a relation- 
ship between n and \A\'^. To this end, one can use the adiabatic equation 
(8). Unfortunately it is impossible to solve this problem analytically for 
the arbitrary temperatures. In the nonrelativistic case (T, Tq ^ 1) Eq. (8) 
gives n = 7(T/To)^/2 and using Eqs. (17)-(18) along with the asymptotic 
expression for G (= 1 -|- |T) for the plasma density we obtain: 



In the ultrarelativistic case (T, Tq :» 1) from Eq. (8) we have n — ^{T/TqY 
and using the asymptotic expression for G (= AT) we get: 




(19) 



The expressions (19) and (20) show that the total density of plasma can be- 
come zero at \A\^ — \A\lj. (where = STq for the nonrelativistic case and 
\A\l^ — 16Tq for the ultra- relativistic case). This phenomenon can be called 
" electron-positron cavitation" and for the ultrarelativistic case has been dis- 
covered in Ref. [8]. It should be noted that the upper bound limitation of the 
amplitude of the vector potential is caused by the fact that we consider here 
the stationary case, and consequently the inertial terms in Eqs.(12)-(13) have 
been neglected. In the case when > \A\l^, a gas-dynamic pressure force 
cannot compensate the ponderomotive one and a stationary distribution of 
fields docs not exist. 

In the ultrarelativistic case, / = |Ap/16rQ and Eq. (15) takes the form 
of the well known nonlinear Schrodinger equation. The stationary soliton 
solution of this equation (which corresponds to nonlinear frequency shift 
A = -A^/16To^) is 

A = Amsech (^z\ (21) 



were A^ is amplitude of soliton. Note that amplitude of the soliton A^ can 
be relativistically strong {Am ^1). The only restriction is that Am < Acr = 
4To. If Am —>■ Acr then the cavitation of plasma occurs and all particles are 
rejected from the central part of soliton. 

Now let us consider the nonrelativistic case. Using Eq. (19) / can be 
written as 

^-'-('-^r '''' 

substituting Eq.(22) into Eq.(15) we get the nonlinear Schrodinger equation 
with a saturating nonlinearity. For the stationary soliton solution we should 



9 



solve to following equation: 

^-\^E + 2E[l - (1 - E'^f'] = (23) 

where E = \A\/(bToYl'^ < 1 and A = -A^. The Eq.(23) has a soliton solution 
provided that nonlinear frequency shift satisfies the following "dispersion" 
relation: 

A^ = 2-^^(l-(l-i^^)^/^) (24) 
where Ej^ is the amplitude of the soliton. One can see that monotonically 
grows with and when E'^ — > 1 (which corresponds to cavitation) obtains 
its maximal allowed value A^^^, — 1.2. Unfortunately the general analytical 
solutions of Eq.(23) can not be expressed in terms of elementary functions. 
We would like to mention that in the case of small amplitude solitons {Em -C 
1) the soliton represents a soliton solution of the cubic nonlinear Schrodinger 
equation and can be written as: 

.3. 1/2 



E = EmSech ^-j EmZ (25) 

If Ejn ^ 1 (A^ — > 1.2) the top part of the soliton is well described by a cosine 
function and can be approximated as 

E ^ EmCos[{2 - XY^h] (26) 

The general shape of the soliton is displayed in Fig.l where E^ = 0.99 
(A^ = 1.19). Dashed line corresponds to analytical approximation given by 
Eq.(26). 

Using the well-known stability criterion of Vakhitov and Kolokolov [9] , it 
can be shown that the above described soliton-type solution is stable against 
small perturbations. 
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In conclusion, we have considered the possibihty of high-frequency EM 
wave locaUzation in hot unmagnetized electron-positron plasmas. In our 
analysis, we included not only the relativistic effects in the hydrodynamic 
motion of the plasma, but also the effects which result from the relativis- 
tic electron velocity distribution. We have shown that in such plasmas it 
is possible to have localized stable soliton-like structures. It is also shown 
that cavitation of plasma can occur both in the nonrelativistic and ultrarel- 
ativistic cases. The present result should be useful for the understanding of 
the nonlinear photon motion in cosmical plasmas such as those found in the 
early universe and AGN. 
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Figure Captions 

Fig.l Solution £^ as a function of the space coordinate {Em — 
The dashed line corresponds to analytical approximation given by Eq. 
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